Transition from the Compact to the Dense Phase of Two-Dimensional Polymers 
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We present a unifying picture of the compact, dense and dilute phases of two-dimensional poly- 
mers. The lattice dependence of the scaling exponents for compact polymers is reconciled with their 
universality in the dense and dilute case. In particular, we show that violations of the fully-packing 
constraint in the compact phase can be interpreted as magnetic screening in the associated Coulomb 
gas, which induces a flow to either the dense or the dilute phase. When more than one flavour of 
polymers is present the flow away from the compact phase leads to a decoupling of the flavours, and 
the central charge decreases by an integer. If charge asymmetry develops the polymer flavours may 
independently flow to either of the two non-compact phases. 



I. INTRODUCTION 



Unlike polymers in the much studied dense and di- 
lute phases |lj , two-dimensional polymers in the compact 
phase exhibit a remarkable lack of universality. Namely, 
their scaling properties are described by conformal field 
theories that depend crucially on the coordination num- 
ber of the underlying lattice. This leads to an unusual 
situation where critical exponents for compact polymers 
on the square and the honeycomb lattice are different 



Another intriguing property of compact polymers, 
which by definition cover all the lattice sites, is that they 
constitute new universality classes that differ from that 
of dense polymers , for which the fraction of covered 
sites < / < 1. In the dilute phase, which corresponds 
to a polymer in a good solvent, / = 0. 

The difference between compact and dense polymers is 
easily understood from the random surface perspective. 
Namely, two-dimensional polymers can be mapped to a 
lattice model of a fluctuating interface which is described 
in the scaling limit by a Liouvillc field theory with imag- 
inary couplings. However, whilst the dense and dilute 
phases are both adequately described by a scalar height 1 , 



the fully-packing (or Hamiltonian) constraint imposed on 
compact polymers forces the height space to be (z — 1)- 
dimensional, z being the coordination number of the lat- 
tice at hand. The larger dimensionality of the height 
space leads to a more complicated effective field theory 
and to new exponents pH. 



Nonetheless, certain critical exponents turn out to be 
"super-universal" , in the sense that they are identical 
for dense polymers and compact polymers on the square 
and honeycomb lattice. Consider the exponents x s , for 
s = 1, 2, . . ., governing the power-law decay of the corre- 
lation function G s (r) ~ r~ 2Xs which describes s polymer 
chains with one end anchored within a small region cen- 
tred around the origin, and the other in a region around a 
distant point r. A comparison of the explicit expressions 
for these so-called string (or watermelon) dimensions for 
the three models under consideration reveals that they 
are the same for s even. This is also true of the differ- 
ence x s — x s > between any two odd string dimensions. 
These observations have important consequences for the 
physically relevant contact exponents 2 Og. These expo- 
nents determine the asymptotic decay of the probability 
Pg{y) ~ y 9 ° as y — > 0, that a certain number of points 
on the polymer separated by macroscopic distances along 
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1 When more than one loop flavour can be defined we shall, 
not surprisingly, find that one height component is needed for 
each of these flavours. However, unlike what is the case for 
compact polymers these height components decouple in the 
effective field theory. 

Globular proteins in their native state form compact struc- 
tures, and hydrophobic interactions that take place at con- 
tacts play an important role in the folding process M . 
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the chain, simultaneously have a spatial separation of or- 
der y; the arrangement of contact points is specified by 
the graph Q. Restricting for simplicity the attention to 
two-point contacts, it is easily seen that 9g only depends 
on the string dimensions through the universal combi- 
nations X3 — x\ and long as only interior points 
and at most one of the endpoints of the polymer are 
involved in the contact. On the other hand, contacts 
involving both endpoints give rise to exponents that de- 
pend on x\ separately, and hence arc not universal. Inci- 
dentally, x\ also determines the conformational exponent 
7 = 1 — x±, which describes the scaling of the ratio of 
the number of open to closed polymer chains with their 
length, and is a measure of the effective entropic interac- 
tion between the two polymer ends. For compact poly- 
mers on the honeycomb lattice the mean-field ]9| result 
xf 1 = was found HQ, whereas both compact poly- 
mers on the square lattice and dense polymers exhibit 
an effective entropic repulsion between chain ends with 
xf = -5/112 and xf = -3/16 respectively §J|. 

This paper aims at providing a theory of the tran- 
sition from the compact to the dense phase, with em- 
phasis on the mechanism that restores the universality 
of the string dimensions when violations of the fully- 
packing constraint are permitted. A particularly inter- 
esting framework for these investigations is an exten- 
sion of the square-lattice model which was used to solve 
the compact polymer problem in Ref. Namely, this 
model allows for the definition of two mutually exclud- 
ing flavours (species) of polymers, which in the compact 
case were found to interact in a manner that endowed 
the model with a two-dimensional manifold of critical 
fixed points. Apart from the obvious interest in exam- 
ining what happens to the Liouville field theory as one 
flows away from the compact phase, it is a most beguil- 
ing question whether the dense and dilute phases also 
allow for a two-flavour extension. We shall find the an- 
swer to this question to be affirmative, with an important 
caveat: Once vertices are allowed to be uncovered by the 
polymers, the two flavours decouple. Quite remarkably, 
this decoupling turns out to be so complete that one of 
the flavours may even flow to the dilute phase whilst the 
other remains dense! This provides a rationalisation for 
four of the five critical branches of the O(n) model stud- 
ied in Ref. 0]. 

We begin (in Section [n]) by giving a heuristic argu- 
ment, based on the field theory for the compact phase, 
that explains how dense exponents may emerge from the 
compact ones when we allow for violations of the / = 1 
constraint. In the Coulomb gas formalism string dimen- 
sions such as the ones defined above are determined by 
inserting a pair of electromagnetic test charges into the 
vacuum, and calculating their interaction energy. Vio- 
lations of the compactness constraint are shown to cor- 
respond to magnetic charges which form a plasma that 
completely screens one of the components of the vector 



magnetic charge of the test particles. As a result the 
compact string dimensions reduce to those of two non- 
interacting flavours of dense polymers. 

Though physically appealing this argument is un- 



sound. We therefore go on (in Section [II) to propose a 



lattice model that exhibits all three phases of polymers — 
compact, dense, and dilute. Exact values of the central 
charge and the critical exponents are then computed from 
a Coulomb gas construction which, although not mathe- 
matically rigorous, is well founded in the renormalisation 
group |l[] . These results neatly confirm the more intuitive 
magnetic plasma picture, and go further in describing the 
interplay between the dense and dilute phase. 

Even though the two-flavoured model on the square 
lattice is our prime concern, the ideas presented here 
are quite general. In particular they apply to the one- 
flavoured model on the honeycomb lattice, and in general 
to loop models defined on arbitrary z-fold coordinated 
regular lattices. This is the subject of the Discussion 
(Section 0), where we also elaborate on the relation be- 
tween our results and those obtained for the O(n) model 
of Ref. M. 



II. COULOMB GAS WITH A MAGNETIC 
PLASMA 

A. Two-flavour fully packed loop model 

In Ref. ^ scaling properties of compact polymers on 
the square lattice were investigated with the help of the 
two-flavour fully packed loop (FPL 2 ) model. This model 
has a two-dimensional region of critical fixed points 
wherein compact polymers are identified with a single 
point; this point also belongs to a line of fixed points 
which is a subset of the critical manifold and it describes 
interacting compact polymers. Since our results for the 
compact-to-dense transition extend to the entire critical 
manifold we shall begin by briefly reviewing the FPL 2 
model. 

First define Q' to be the set of all edge colourings of 
the square lattice with four different colours. In other 
words, each edge of the lattice is assigned a colour A, B, 
C or D, subject to the constraint that for any vertex all 
of its four adjacant edges carry different colours. If we 
bipartition the square lattice into even and odd sublat- 
tices, directed loops of two different flavours ("black" and 
"grey") can be defined as follows: Colour A designates 
a black loop segment directed from an even to an odd 
site, and colour B corresponds to a black loop segment 
having the opposite direction. Similarly colours C and 
D are used to define directed grey loops. By Q we denote 
the equivalence classes of Q' with respect to independent 
changes of directions for each of the loops. The FPL 2 
partition function is then 
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z = H< b <^ CO 

Q 

where rib,n g are the fugacities and Nb, N g the numbers 
of (undirected) black and grey loops in Q respectively. 
The model ([[]) is critical on the manifold < rib, n g < 2, 
and the compact polymer problem discussed in the In- 
troduction is recovered in the limit (rib, n g ) — > (0, 1). 

The extra orientational information present in Q' al- 
lows for a local redistribution of the loop fugacities. In- 
deed, assign a phase factor exp(±i7reb/4) to each vertex 
where a directed black loop makes a right (left) turn, 
and the factor 1 if it continues straight. In this way 
the weight of the entire black loop becomes exp(±i7reb) 
if the orientation is clockwise (counter-clockwise). The 
undirected loop weight rib is obtained by summing over 
the two orientations, whence rib = 2cos(7reb). With a 
similar convention for the grey loops, and letting A(x) 
be the product of the local black and grey weights at 
the vertex x, the partition function, Eq. (|l|) , can now be 
rewritten as 

Z = J2U X ^- ( 2 ) 

G' x 

The continuum limit of the FPL 2 model is obtained 
by mapping the oriented loop configurations to an in- 
terface model defined on the dual lattice. The micro- 
scopic heights z live on the plaquettes of the square lat- 
tice and they are defined, up to an overall shift, by the 
edge colours which stipulate the height differences be- 
tween neighbouring plaquettes. More precisely, when en- 
circling an even (odd) site in the clockwise direction a 
vector A, B, C or D is added to (subtracted from) z 
depending on the state of the edge being crossed. The 
fully-packing constraint imposes the condition 

A + B + C + D = 0, (3) 

whence the colour vectors are three-dimensional. A suit- 
able representation is 

A=(-l,+l,+l), B = (+l,+l,-l), 

C = (-1,-1,-1) , D = (+l,-l,+l) . (4) 

The continuum limit is then taken by coarse graining z. 
The effective field theory is described by a Liouville ac- 
tion containing three terms : 

S = Se + Sb + Sb . (5) 

The elastic part <Se is a gradient-squared term in the 
coarse grained height h(x), and it controls the fluctua- 
tions of the interface away from flat configurations. The 
boundary term Sb assigns the correct weights to loops 
winding around the point at infinity by coupling h to the 
scalar curvature 1Z. Finally, Sl is a Liouville potential 
that assigns the correct weights to loops in the bulk, and 



it is the coarse grained version of the microscopic vertex 
weights A(x). 

In the dual Coulomb gas picture the action, Eq. (||) , is 
interpreted as that of a two-dimensional gas of electro- 
magnetic vector charges interacting via a logarithmic po- 
tential. Electric charges are associated with vertex op- 
erators exp(ie • h) which appear as the scaling limits of 
FPL 2 operators that can be expressed as local, periodic 
functions of the microscopic heights. Magnetic charges 
m act as non-local constraints on the height field, and 
they are associated with topological defects in the mi- 
croscopic heights or, equivalently, with violations of the 
four-colouring constraint. Finally, the curvature term 
Sb introduces a background electric charge — 2e at the 
boundary of the system. For reasons of neutrality this 
endows the Coulomb gas vacuum with a floating charge 
+2eo located in the bulk of the system. 

B. Height defects and magnetic screening 

Following the standard Coulomb gas recipe Q , the cal- 
culation of the scaling dimension of an electromagnetic 
operator reduces to the determination of the interaction 
energy of "test particles" of charge ±(e,m) inserted into 
the vacuum. This energy varies logarithmically with sep- 
aration, where the prefactor of the logarithm is twice the 
sought dimension. As the floating charge is free to move 
around the bulk it will coalesce with one of the test par- 
ticles thus minimising the energy. Therefore the interac- 
tion energy consists of three terms associated with the 
pairwise interactions of the background electric charge 
(— 2e ,0), (e, m), and (2e — e, — m), with the vacuum 
energy subtracted off. The energy of the vacuum is sim- 
ply due to the pair interaction between the background 
charge and the floating charge. 

Consider, for example, the computation of the string 
dimensions x s for the FPL 2 model. Since the defect 
strings may now have any of the two flavour labels this 
generalises to x Sb ,s g , where Sb (s g ) denotes the number 
of black (grey) strings. For dislocations generated in the 
bulk the number Sb+Sg must be even Q . Appropriate de- 
fect configurations may be specified through the colour- 
ing state of the four edges around some fixed vertex [jl2| . 
For instance, (A, A, C,D) generates two directed black 
strings and corresponds to a height vortex of strength 
m 2! o = A B = (-2, 0, 2). Similarly, (A, B, C, C) is a 
vortex of strength mo^ = C — D = (—2, 0, —2) that gen- 
erates two grey strings, and (A, C,C,D) generates one 
string of either flavour and has strength rxii.i = C — B = 
(—2,-2,0). For simplicity we shall concentrate on the 
black strings and define x s = x SjSg , where s g = s mod 2. 
The defect magnetic charge 

_ J (— s, 0, s) for s even, , . 

ms ~ \ {-s - 1, -2, s - 1) for s odd, [b) 
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found as a suitable linear combination of ni2.o and mil, 
must be accompanied by an appropriate compensating 
electric charge that takes care of the spurious phase fac- 
tors arising from the winding of the strings around the 
defect cores [El. Setting the fugacity of the grey strings 
equal to unity (e g = 1/3) we obtain the following result 
for the scaling dimensions: 



2x c 



|(1 -e b > 2 
|(l-e b ) S 2 



l-e b 

eg 
l-e b 



for s even, 
2 ^~ eb) for s odd. 

5— 3eb 



(7) 



As mentioned in the Introduction the even scaling dimen- 
sions agree with those of dense polymers [pj (and also 
with those of the fully-packed loop model on the honey- 
comb lattice [Q), whereas the odd dimensions contain an 
extra term due to the coupling between the black and 
the grey strings. 

So far our discussion has been for the fully-packed case 
where a fraction / = 1 of the lattice vertices are visited 
by both loop flavours. Now consider moving infinitcsi- 
mally into the dense phase by diminishing the fraction of 
visited vertices to / = 1 — e. One can imagine doing so 
by distributing a fraction e of (A, B, A, B) defects and a 
fraction e of (C, D, C, D) defects throughout the system. 
The corresponding vortices have magnetic charge ±m<r, 
where 



m T = (0,4,0) 



(8) 



Physically this situation corresponds to having a dilute 
magnetic plasma superimposed on the Coulomb gas vac- 
uum. 

As before we can then imagine inserting two electro- 
magnetic test particles of charge ±(e, m) into the system 
in order to measure the string scaling dimensions. Just as 
in the familiar case of electric test particles in a dielectric 
medium there will now be some screening, the difference 
being that here it is the magnetic charge that is being 
screened. However, unlike what is the case in this dielec- 
tric analogy, the magnetic plasma is completely free to 
move around, and the screening will be complete. This 
is a consequence of the fact that the magnetic charges in 
Eq. (pf) are relevant, i.e. their scaling dimension |q] 



2x(±m T ) = 4 



(l-e b )(l-e 8 ) 
(1 - e b ) + (1 - e g ) 



(9) 



is less than 2 on the whole critical manifold of the 
FPL 2 model. Therefore these charges are not bound into 
vortex-antivortex pairs at large scales Jl3| . 

Taking our cue from this intuitive picture we hypoth- 
esise that the effect of the magnetic plasma can be mod- 
eled by the substitution 



m : 



(to , to 2 , to 3 ) 



m 



(to\0,to 3 ), (10) 



since the only non- vanishing component of the screening 
charge ([|) is along the 2-direction, and this does not cou- 
ple to any of the other directions. In particular, replacing 



m s of Eq. (||) by m' s we find that instead of Eq. (Q) we 
get 



2x' 



1 



(1 - e h )i 



1 - e b 



(11) 



for all s, whether even or odd! It seems that the interac- 
tion between the two loop flavours that was responsible 
for the difference between the even and odd string di- 
mensions in Eq. (^) has somehow been disposed of, and 
what remains is the critical exponents of dense polymers 

One can proceed to compute the central charge for 
the screened system. We recall that in the compact 
(/ = 1) case this was given as c = 3 + 12a;(e , 0), where 
the three free bosons present in Se each contribute unit 
central charge, and there is a negative shift due to the 
background electric charge. However, since the magnetic 
charges live in height space the complete screening of 
to 2 means that height fluctuations along the 2-direction 
are frozen out. This is nothing but the usual Kosterlitz- 
Thouless scenario 13 , where due to the unbinding of 
dislocations the interface defined by the second height 
component is rendered smooth. As a result the central 
charge for / = 1 — e is d = c— 1 = 2 + 12a;(eo, 0'), which 
can be written as 



c'(e b ,e g ) = C dn ( eb )+c d "(e g ) , 



where 



„dn, 



e) = l 



6e 2 
1 - 



(12) 



(13) 



is nothing but the central charge of one flavour of dense 
polymers 0. If we again set the fugacity of the grey 
loops equal to unity (e g = 1/3) we find that c dn (e g ) = 0, 
and only the contribution from the black loops remains. 
The sum rule (12) should be taken as a further indica- 
tion that once we allow for violations of the fully-packing 
constraint the two loop flavours decouple. 

At this point, of course, a number of objections may be 
raised. First, it is seen from Eq. (^|) that the magnitude 
of the magnetic charge along the 2-direction is —2. How- 
ever, for the screening to take place it then seems that 
the magnetic screening charges im^ given by Eq. (||) 
must somehow be fractionalised. 

Another objection is that the above calculation, when 
taken at face value, purports to be perturbative around 
the fully-packed phase. However, as pointed out earlier 
[Eq. M)] the defect operators within the FPL 2 model are 
strongly relevant! This means that the fugacity of these 
defects will flow off to infinity at large scales. Hence, 
the magnetic plasma analogy is no more than a heuristic 
argument. Nevertheless, we shall see in the following sec- 
tion that it is possible to generalise the FPL 2 model to 
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a lattice model that explicitly accommodates configura- 
tions that are not fully packed. In this generalised model 
Coulomb gas calculations can be carried out, confirming 
the main results found in this Section. 

Even when accepting the magnetic plasma analogy as 
a heuristic framework for treating the case of / < 1, 
the FPL 2 height mapping on which it is based breaks 
down completely for / — > 0. It is therefore impossible to 
say anything about the role of the dilute polymer phase. 
However, from the calculations presented in the next Sec- 
tion we shall be able to say more about how dilute expo- 
nents may emerge. 



III. TWO-FLAVOUR DENSELY PACKED LOOP 
MODEL 



In order to unify our understanding of compact, dense 
and dilute polymers as well as their possible interrelation 
we introduce a statistical mechanics model that can ac- 
commodate all three phases. Since the critical properties 
of the compact phase are known to be lattice dependent 
we need to consider the model on a specific lattice. To ex- 
amine the important issues associated with the presence 
of more than one loop flavour we choose the square lat- 
tice, whereas the question of what the construction would 
look like on other lattices is deferred to the Discussion. 



h r 4 4 



(14) 



where rib, %j -Wb> N g have the same meaning as in Eq. (Q), 
and Vb and v g are the number of occurrences of vertices 
7-8 and 9-10, respectively. The summation runs over 
the set Q A of loop configurations on the square lattice 
that can be made out of vertices 1-10, subject to the 
constraint that both flavours form closed loops. This 
model is closely related to the O(n) model introduced in 
Ref. [|l0| for which five branches of critical points were 
found as a function of n. In the Discussion we shall com- 
ment on the precise relation of the DPL 2 model to this 
one, thus providing a physical picture of its critical prop- 
erties. 

Clearly, in the special case (W h , W s ) = (0, 0) the DPL 2 
model reduces to the FPL 2 model, and the analysis of 
Ref. H applies. Let us briefly recall the main results of 
this fully packed case. For Sb + s g = mod 2 the full 
spectrum of string dimensions is 



2z Sb , Sg = \ [(l-e b )s 2 + (l-e g )s 2 ] 



1 - e b 



1 



(15) 



r(2) c(2) (l-e b )(l-e g ) 

^V (1 _ et) + (1 _ eg) : 



where 6\ 2 } 
charge is 



Ji=j (mod 2) ■ 



c = 3-6 



1 — eb 



The result for the central 



(16) 
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8 



10 



FIG. 1. The ten allowed vertices in the DPL 2 model. Black 
and grey loop segments are represented by solid and dashed 
linestyle respectively. 

We therefore define the two-flavour densely packed loop 
(DPL 2 ) model on the square lattice by listing its ten al- 
lowed vertex states in Fig. |l|. The first six vertices are 
those familiar from the fully packed (FPL 2 ) case, whilst 
the last four vertices explicitly allow the model to vio- 
late the fully-packing constraint. Vertices 7-8 exclude 
the grey loops from a given site and carry a weight Wb, 
and similarly vertices 9-10 exclude the black loops and 
have weight W g . The first six vertices are all assigned 
unit weight. The partition function of the DPL 2 model 
is then 



In the following we focus on the Coulomb gas construc- 
tion for the case where either Wb or W g is greater than 
zero. 



A. Colouring and height representation 

The mapping from oriented loop configurations Q' A to 
the four-colouring representation goes through exactly as 
described for the FPL 2 model in Sec. II A. When defining 
the microscopic heights z there is however one very im- 
portant difference. The presence of vertices 7-10 in Fig. |l| 
imposes extra constraints on the choice of the colour vec- 
tors: 



A + B = 



C + D = 



(17) 



Actually, whenever one of these constraints holds true 
the other follows from Eq. ([|). So the conditions ( p^ ) 
are valid provided that either Wb > or W g > 0. 

As a consequence the color vectors in Eq. (^) must be 
replaced by 
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A = (+1,0), B = (-1,0), 
C = (0,+1), D = (0,-1), 



(18) 



i.e., the height space is now two-dimensional. This elim- 
ination of one of the height components is very reminis- 
cent of the effect of screening in the magnetic plasma 
analogy of Sec. O. 

The local redistribution of the loop fugacities, and 
n g of Eq. (|l4|) , in terms of local weight factors associated 
with the left an d rig ht turns of the loops, again takes 
place as in Sec. II A. One may wonder how should the 
vertex weights Wb and W g be represented? Unlike the 
loop weights which do not renormalise under a change of 
scale 3 , the vertex weights do change under scale trans- 
formations as they flow towards their fixed-point (and 
non-universal) values. Thus, from the point of view of 
the Liouville field theory description that we are about to 
set up for the scaling limit, they need not be represented 
at all 4 ! Of course, all along our working assumption is 
that the scaling limit of the DPL 2 model exists for loop 
fugacities < rih,, n g < 2 . 



B. Ideal states 



An ideal state of a loop model has two distinguish- 
ing features. First, it is macroscopically flat, so that 
fluctuations of the interface around the ideal states can 
be controlled by an elastic term 5e in the effective ac- 
tion. Second, it is entropically selected in the sense 
that it maximises the number of configurations which 
differ from it by the smallest allowed change. To be 
concrete consider the directed FPL 2 model in the four- 
colouring representation. Since any change of a con- 
figuration must necessarily involve the change of some 
edge colour o\ into another colour 02, it follows that 
the whole loop o\aio~\o~2 ■ ■ ■ o\Oi defined as an alternat- 
ing sequence of these two colours must be changed into 
02<J\<Ji<J\ ■••<72<Tx. An ideal state allows for the maxi- 
mum number of such loop flips, and consequently has a 
maximum number of small loops of alternating colour. 



I 



(a) (b) (c) 2 

FIG. 2. Representative examples of ideal states in the FPL 
model, here shown in the oriented loop representation. The 
full ideal state is obtained by tiling IR 2 with the 2x2 pattern 
given in the figure. Each of the examples corresponds to a 
group of eight ideal states, related by independent changes of 
orientations for the two flavours, and by flavour conjugation. 



In the FPL 2 model on the square lattice there are 24 
ideal states, defined by the 4! possible permutations of 
the colours A, B, C and D around a fixed vertex. For 
later convenience we choose to label them by the clock- 
wise sequence (o\, 02, 03, 04) of colours around this ver- 
tex, starting with the leftmost edge. In the oriented loop 
representation the ideal states form a pattern of 2 x 2 pla- 
quettes which is repeated throughout the lattice. Three 
representative examples are shown in Fig. ||.a~c. Taking 
into account the possible changes of loop flavours and di- 
rections there are eight FPL 2 ideal states corresponding 
to each of these examples, making a total of 24. 



When we allow for the extra vertices of the DPL 2 
model more ideal states become possible. Namely, for 
each of the 16 states exemplified by Fig. ||.a-b there cor- 
responds another ideal state where only one of the two 
loop flavours is used to form the small loops. Again we 
can label the states by means of the colour configuration 
around a fixed vertex, but as shown in Fig. |^ this is in 
some cases no longer unique. To remedy this we adopt 
the convention to add a prime to the label whenever the 
loop configuration around a vertex is as shown in Fig. |[a; 
see Fig. |[ 



3 This is a consequence of the loop ansatz pj],|q| to which we 
shall return in Section 
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below. 



4 Actually, this is not quite true. As mentioned above the 
constraints (^) explicitly take into account that at least one 
of Wh or W s is non-zero. Furthermore, one may suspect that 
if the bare value of Wh is very large the grey loops flow to 
the dilute phase. This is indeed the case, and we shall see 
below how this mechanism is implemented in the Coulomb 
gas construction. 
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(B,A,B,A) BAB A' BABA 

FIG. 3. Certain colour configurations around a vertex do 
not uniquely specify the directed loop configuration. In par- 
ticular the label (B,A,B,A) corresponds to two distinct 
DPL 2 ideal states, which we distinguish by a prime as shown. 

This ambiguity has also bearings on the way in which 
changes of a given colouring configuration may be per- 
formed. Namely, when changing some edge colour o\ 
into (72 one may encounter one of the ambiguous ver- 
tices of Fig. U in the process of tracing out the loop 
0i 02 <7 1C2 • • • 01(72 to be flipped. In that case one has 
two equivalent choices for proceding. However, this pe- 
culiarity does not change the fact that the 16 new states 
are still ideal, since they are perfectly macroscopically 
flat and they allow for an equal number of small changes 
as do the 24 ideal states of the FPL 2 model. The DPL 2 
model thus has 40 ideal states. 

The lesson to be learned from this complication is that 
it is favourable to think of the DPL 2 ideal states in terms 
of the oriented loop representation rather than in terms 
of the colouring representation. In fact, the mapping 
from oriented loop configurations to the four-colouring 
representation is now many-to-one, so that the inverse 
mapping is no longer well-defined. Since we only need 
the colours to define the heights, which are the basic con- 
stituents of the continuum field theory, we do not have 
to worry whether the series of mappings can be reversed. 



C. More on the dimensionality of height space 



It is worthwhile noticing that the choice ( |l8| ) would 
have been just as good as (|]) in ensuring the one-to- 
one correspondence between oriented FPL 2 configura- 
tions and the microscopic heights. However, a too re- 
strictive choice of the colour vectors may render a non- 
ideal state macroscopically flat, in which case the elastic 
term 5e of the action fails to enforce the proper entropic 
penalty for a state that does not allow for a maximum 
number of loop flips. 
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FIG. 4. This FPL 2 state has a non-vanishing macroscopic 
height gradient along the horizontal direction. It is shown in 
three equivalent representations: (a) oriented loops, (b) edge 
colourings, and (c) microscopic heights. 

As an example consider the state arising from a tiling of 
IR 2 with the 2x8 pattern shown in Fig. [|. Upon coarse 
graining this state is flat in the vertical (x 2 ) direction, 
but has a macroscopic height gradient in the horizontal 
direction: 



1 



dh/dx 1 = - [(C + D) - (A + B)] 



(19) 



With the choice ( |l8| ) the height gradient vanishes, and 
5e = for this height configuration. On the other hand, 
the choice @ leads to dh/dx 1 = (0,-1/2,0), and the 
state is duly suppressed by Se since the stiffness con- 
stant along the second height direction 5 



K 22 = 



(l-e b )(l-e 8 ) 
(1 - e b ) + (1 - e g ) 



(20) 



is non-zero throughout the critical region 

This important point is further illuminated by noting 
that black and grey loops play the role of contour lines for 
the projection of h along the (1,0, —1) and the (1,0, 1) 
directions respectively. Evidently these projections fail 
to detect any variation along the 2-direction in height 
space, and accordingly the state given by Fig. |^ appears 
macroscopically flat when viewed exclusively in terms of 
its oriented loop representation. 



5 The similarity between this expression and Eq. (^) for the 
thermal scaling dimension emphasises the special role of the 
second height component. 
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D. Ideal state graph 



the truncated octahedra. In particular, the square and 
the hexagonal faces are clearly visible. 



The 40 ideal states of the DPL 2 model define the ideal 
state graph X which is instrumental in the construction 
of the continuum field theory. The idea is to identify the 
ideal states with nodes on X whose position is given by the 
average microscopic height (z) of the given state. Two 
ideal states are said to be nearest neighbours if their re- 
spective microscopic heights are identical on three fourths 
of the lattice plaquettes. On X, the two nodes corre- 
sponding to such a nearest neighbour pair are connected 
through a link. 




FIG. 5. Ideal state graph 1 of the DPL model. The thick 
and thin linestyles distinguish the links which are inherited 
from the FPL 2 model from those which are new to the DPL 2 
model. The scale is set by the squares in thick linestyle, which 
have side length 1/2. 

The result of this construction is the ideal state graph 
shown in Fig. |[ Here the links added due to the inclusion 
of the 16 ideal states which are new to the DPL 2 model 
are shown as thin lines, whereas the thick linestyle desig- 
nates links between the 24 ideal states which have been 
inherited from the FPL 2 model. We recall that the ideal 
state graph for the FPL 2 model was a covering of IR 3 by 
truncated octahedra. Since the DPL 2 colour vectors ( p"§| ) 
can be obtained from their FPL 2 counterparts (Q) by a 
projection onto the z 1 z 3 -plane and a subsequent ^ ro- 
tation and rescaling, it should hardly come as a surprise 
that the thick lines form a two-dimensional projection of 




FIG. 6. The repeating unit of 2" along with the state la- 
bels. Ideal states that are new to the DPL 2 model can easily 
be distinguished by noting that only two different colour vec- 
tors appear in their state label. 

As a consequence of our definition of nearest neighbour 
ideal states any one ideal state corresponds to an infinity 
of nodes in X. The entire ideal state graph thus corre- 
sponds to a covering of IR 2 by the repeating unit shown 
in Fig. |(| Two novel features of the DPL 2 ideal state 
graph which have not been encountered in the previously 
studied loop models are that each node may have several 
labels, and that links may cross without being connected 
by a node. The latter observation just indicates that X 
really must be thought of as a graph (a collection of nodes 
and links) and not simply a lattice; of course, unlike the 
usual notion of a graph, the precise location of the nodes 
is of paramount importance as it bears on the effective 
field theory of the model. 

Concentrating first on the 24 ideal states on Fig. |^ that 
are connected by thick links the division of the states 
into the three classes exemplified by Fig. ^| becomes ob- 
vious. The two groups of eight states having the ele- 
mentary plaquette loops positioned as on Fig. ^|.a and 
Fig. |^.b respectively, each form a square with doubly la- 
beled nodes. These squares are connected by thick lines 
via two four- fold labeled nodes ( "crosses" on Fig. |) that 
together accommodate the eight states of the type shown 
m Fig. |.c. 6 The remaining 16 ideal states which are new 



3 At first sight these states hardly appear "ideal", at least 
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to the DPL 2 model, located at intersections of thin lines 
in Fig. occur in two groups of eight states correspond- 
ing to plaquette loops positioned as in Fig. |^.a (primed 
states) and in Fig. |^.b (unprimed states). Quite natu- 
rally, the primed states are linked to the square of eight 
FPL 2 states with loops positioned as shown in Fig. |[a, 
and similarly the unprimed states are distributed around 
the other square (c/r. Fig. ||.b). In particular, the states 
DCDC and DCDC' must be thought of as labeling two 
distinct nodes, which just happen to occupy identical po- 
sitions on Fig. |6|. 7 This subtlety becomes more evident 
if one tries to trace out X explicitly by moving from one 
ideal state to the neighbouring one. This is accomplished 
by either choosing two colours in the starting ideal state 
and then interchanging them along all the alternatingly 
coloured loops to obtain a new state, or by changing one 
loop flavour to the other (this leads to a new ideal state 
only for states that are made up of all gray or all black 
loops). It then turns out to be impossible to move be- 
tween states on two different squares using only the links 
that are shown as thin lines on Fig. 0. 



E. Liouville field theory 

A continuum field theory for the DPL 2 model can now 
be obtained by coarse graining the microscopic heights 
over domains of ideal states [O] . When defining the con- 
tinuum height field, nodes in X representing the same 
ideal state should be identified, which means that the 
height must be compactified with respect to the repeat 
lattice 1Z, 



h(x) e \R 2 /1Z . 



(21) 



According to Fig. |^ the repeat lattice is a square lattice 
of side \/2, spanned by the vectors (1,1) and (1,-1). 

We are now ready to write down the Liouville field 
theory for the height field h(x). As usual the partition 
function Z > , which describes only the large-scale fluctua- 
tions of the height, can be written as a functional integral 



Z > = Jvh exp(-5[h]), 
S = Se + Sb + Sl, 



(22) 



where the effective Euclidean action consists of three 



terms that were mentioned in section [I A . Here we dis- 
cuss them in some detail. 



1. Elastic term 



The most general form of the elastic term has the form 



(23) 



where Latin letters label the basal plane coordinates and 
Greek letters pertain to the height space. Since, by defi- 
nition of the allowed DPL 2 vertices (see Fig. [l]), the con- 
tinuum model is invariant with respect to rotations in the 
basal plane, the stiffness tensor is diagonal in the Latin 
indices: K^g = K a pS lJ . The remaining four compo- 
nents are constrained by colour symmetries, correspond- 
ing to independent reversals of the two flavours of di- 
rected loops: 



A <-> B 
C <-> D 



-eb and z 1 
-e„ and z 2 



(24) 



Any one of these symmetries prevents the term K^dh 1 ■ 
dti 2 from occurring in Se- The elastic term thus assumes 
the diagonal form 

S E = \ J d 2 x [K^dh 1 ) 2 + K 2 (8h 2 ) 2 } , (25) 

involving only two independent stiffness constants K\ 
and K-2,. 



2. Boundary term 
The boundary term Sb generically takes on the form 

5b = i / d2x ( e ° • ' ( 26 ) 

where 1Z is the scalar curvature, and eo is the background 
electric charge which is to be determined. To calculate 
eo it is most convenient to consider the DPL 2 model de- 
fined on a cylinder, so that the curvature is concentrated 
at either end: 



U = 4tt[£(+oo) - S(-oo)] 



(27) 



This terms in the action has the effect of placing two 
vertex operators with charges ±eo at the far ends of the 
cylinder. They in turn give complex weights to loops 



not when viewed solely in the oriented loop representation. 
However, it should be clear that they are needed to ensure 
the connectivity of X, c/r. Fig. |H[ 

7 A similar remark holds true, of course, for the primed and 
the unprimed versions of states ABAB, BABA and CDCD. 
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winding around the cylinder. In order for these to be the 
same as those defined earlier for loops in the bulk, the 
equations 



e • A = 7re b , 
e • C = 7re g , 



-7re b 



e B 



ep • D = — 7re 



g ' 



must be satisfied. The unique solution is 
e = (7re b ,7re g ) . 

3. Liouville potential 



According to Eqs. (||) and 
of the Liouville potential is 



(28) 
(29) 



2J) the appropriate form 



St, = 



;[h(x)] 



(30) 



where exp(— w[h(x)]) is the scaling limit of the vertex 
weights A(x) arising from the local redistribution of the 
loop fugacities n b and n g . Since these are uniform in each 
of the ideal states we can consider w(h) to be a function 
of the coarse grained height hel. 

It is reassuring to check that whenever a node in 
X carries multiple labels (cfr. Fig. ||) all of the con- 
cerned states have the same value of the vertex weight 
A(x). The reason is that the corresponding oriented loop 
configurations either (a) are connected by a translation 
in the basal plane, or (b) have a local cancellation of 
left and right turns, or (c) all carry unit weight. The 
only apparent exception is pairs of states of the type 
DCDC, DCDC'. These pairs have already been dis- 



cussed in Sec. [II D, where we found that the two states 



should really be associated with two distinct nodes, which 
just happen to be spatially coincident. The fact that 
w(DCDC) = +me g /2, whereas id(DCDC') = -iwe g /2 
therefore does not lead to any contradiction. 

Any local function of the colours (operator) which is 
uniform in the ideal states defines a periodic function of 
h G T with the periods forming the repeat lattice 1Z. 
Upon coarse graining of the height this periodicity prop- 
erty is conserved. Therefore the continuum limit of any 
such operator can be expanded as a Fourier series, where 
the basis functions in the expansion are the vertex oper- 
ators exp(ie ■ h(x)). Consequently, 



w[h(x)] = 2_, w e exp(ie • h(x)) 



(31) 



e£R* 



Note that the sum does not run over the reciprocal lattice 
1Z* but rather over some sublattice TZ^ C 1Z* . This is be- 
cause w(h) may have some additional periodicity within 
the unit cell of I shown in Fig. ||. A careful analysis re- 
veals that this is indeed the case: The height periods of 
w(h) form a square lattice spanned by (1,0) and (0,1), 
and so TZ^ is a square lattice of side 2-7T, spanned by 
(2tt,0) and (0,2tt). 



4- Loop ansatz 



In order to relate the stiffness constants K\ and K-i 
to e b and e g we need only determine the most relevant 
vertex operators in the expansion ( |3l| ) and impose the 
condition that they be exactly marginal. This is the 
loop ansatz which ensures that loop weights do not renor- 
malise under a change of scale [ fPH . From the viewpoint 
of conformal field theory the electric charges appearing 
in these most relevant vertex operators play the role of 
the screening charges introduced by Dotsenko and Fateev 
©• 

The scaling dimension of a general electromagnetic op- 
erator is given by fl~5|,|6| 



2x(e,m) = — |e • K" 1 • (e - 2e ) +m • K • ml , (32) 



where the stiffness tensor K according to Sec. Ill E 1 as- 
sumes the form 



Ki 
K 2 



(33) 



The candidates for the screening charges are the four 
shortest vectors in 1Z* 



3« = (+2tt,0), e£> = (-27r,0), 
i 3) = (0,+2tt) , e&> = (0,-270, 



(34) 



and the scaling dimensions of the corresponding vertex 
operators are 



z(ei 13 ,0) = 



tt(1 - e b ) 
tt(1 - e b ) 



*(el 2) ,0) 



7r(l + e b ) 



Kel 4 \0) = ^±^. (35) 



Ko 



At (n b ,n g ) = (2,2) the background electric charge van- 
ishes and the four operators are pairwise degenerate. 
However, for a generic point on the critical manifold 
< n b ,n g < 2 the operators exp(ieiP • h(x)) and 

(3) 

exp(ieiu ■ h(x)) are more relevant than the other two. 
The loop ansatz now leads to the result 



tfi = f(l-e b ), K 2 = ^(l-e s ) 



(36) 



In Sec. IIIH we shall see that this choice for the elas- 



tic constants reproduces the critical exponents for two 
non-interacting flavours of dense polymers. 
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F. Charge asymmetry and the dilute phase 

In the case of a Coulomb gas with scalar electric 
charges it has been shown by Nienhuis 0] that for partic- 
ular values of the bare particle fugacities a charge asym- 
metry may evolve, so that the renormalised Coulomb gas 
contains unit charges of one sign only. Unit charges of 
the opposite sign then appear only as excitations, having 
as a result that the amplitude of the most relevant term 
in their two-point correlator vanishes. 

An inspection of Nienhuis' argument reveals that this 
mechanism is rather general, and in particular it may 
occur in the case of electromagnetic vector charges that 
we are concerned with here. In our Liouville field theory 
approach the vanishing of the above-mentioned ampli- 
tude translates into the possibility of havin g on e or more 
vanishing expansion coefficients w e in Eq. (31). 

Obviously it would be quite difficult to determine for 
which values of the bare vertex weights Wb and W g such 
a charge asymmetry may evolve. Here we shall just men- 
tion that the four candidates for the screening charges 
given in Eq. ( |34| ) in general warrant the following possi- 
bilities for the clastic constants 



K l = -{lTe h ) 



K> 



-(lTe g ) . 



(37) 



As already mentioned the upper signs will lead to the crit- 
ical exponents of dense polymers, and they correspond to 
the charge symmetric case. A remarkable observation is 
that the lower signs similarly reproduce the critical ex- 
ponents of two decoupled flavours of dilute polymers 8 . 
Furthermore, nothing seems to prevent a charge asym- 
metry from developing which would select the upper sign 
for one of the couplings and the lower one for the other. 
Thus, the complete decoupling of the loop flavours allows 
for a situation where they each reside in either of the two 
non-compact phases {i.e., dense or dilute). 



G. Symmetry algebra 

The FPL 2 model at the point (rib,"g) = (2,2) was 
shown to possess, in the continuum limit, an su{A) affinc 
Lie algebra symmetry, at level k = 1 [EiJ. This follows 
from the effective field theory of the model which is given 
by a three component height (three free massless bosons) 
compactified on the root lattice of su(4). (As in any con- 
formal field theory there are actually two copies of the 
symmetry algebra, associated with the holomorphic and 
the antiholomorphic components of the height field.) Yet 
another view of what transpires when the vertices that 



allow for the violation of the fully packing constraint are 
included, is provided by the question: What happens to 
the symmetry algebra? 

In order to address this question we examine the DPL 2 
model at the point (n^, n g ) — (2, 2). In previous sections 
it was shown that the effective field theory of this model 
is given by a two-component height with the two com- 
ponents completely decoupled. From the analysis of the 
ideal state graph it was further shown that each com- 
ponent is compactified on the one dimensional lattice 
1Z = \/27.. The combination of the lattice constant and 
the calculated stiffness K = K x = K 2 = tt/2 (Eq. @)) 
is such that this lattice can be identified with the root 
lattice of an su(2) algebra. Therefore, the DPL 2 model 
at the special point is characterised in the scaling limit 
by an su(2)k=i S>su(2)k=i affine Lie algebra. The heights 
make up the free field representation of the appropriate 
Wess-Zumino-Witten model. 

The generators of the holomorphic half of the symme- 
try algebra are the modes of the holomorphic currents 
p7| . There is a total of six currents, i.e., conformal di- 
mension (1, 0) operators, three for each of the two su(2) 
algebras: 



\dh\z), exp(ieW-h), exp(iei 2 » ■ h(z)) 



(38) 



and 



idh^(z), exp(ie^ • h(z)), exp(ie^ • h(z)) , (39) 

where the electric charges appearing in the above vertex 
operators are the screening charges of Eq. (pih. Here 



z = x + ix is the complex coordinate, whilst h(z) 
is the holomorphic component of the height field, i.e., 



h(a; 1 



h(z) + h(z). 



Finally we can give a concise answer to the question 
posed at the beginning of this Section: The symmetry 
algebra of the compact phase is lowered from su(4) to 
su(2) © su(2), as it flows to the dense phase. One of 
the three free bosons of sw(4) becomes massive, whilst 
the other two each form a free field representation of an 
su(2)fc = i affine Lie algebra. 



H. Critical exponents 

Having expressed the elastic constants in terms of the 
loop fugacities, Eq. (|37|), we are now in the position to 
compute the central charge and the critical exponents of 
the DPL 2 model in the critical region < nt,,n g < 2. 



8 A similar phenomenon is observed in the Liouville field the- 
ory solution of the O(n) model flfl]. 
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1. Central charge 



IV. DISCUSSION 



For (rib,n g ) = (2,2) the background charge vanishes, 
and the Liouville theory (|2^) consists only of the elas- 
tic term Sr. This is nothing but the action of two free 
bosons, yielding a central charge of 2. In the general case 
there will be a shift due to the background electric charge 




c = 2 + 12.x(e ,O) 
= 2 - 6 f T ^- + -^-j. (40) 

If one chooses the upper signs this result agrees with 
that of the magnetic plasma analogy, and according to 
Eq. (|l2|) it is consistent with a situation where the two 
loop flavours have decoupled and are both in the dense 
phase of the O(n) loop model. More generally we can 
write 

c(e h ,e g )=c dn - d \e h ) + c dn - d \e g ) , (41) 

where c dn ' dl (e) = 1 — 6e 2 / (1 =p e) is the central charge of 
a single flavour of dense (dilute) loops. 



In this paper we have presented a general mechanism 
by which compact polymers may flow to the dense and di- 
lute phases once violations of the fully-packing constraint 
are allowed. The universal dense exponents were re- 
trieved from the lattice-dependent compact ones through 
a heuristic analogy to a magnetic plasma, or alterna- 
tively by a Coulomb gas solution of a statistical mechan- 
ics model (the DPL 2 model) where two adjustable pa- 
rameters (Wb an d W g ) were used to control the devia- 
tions from fully packing. We have found that the com- 
pact phase is unstable to any such deviation, however 
small, and that the flow may be to the dense or the di- 
lute phase depending on the bare values of the adjustable 
parameters. Furthermore, by this mechanism the two 
loop flavours of the FPL 2 model decouple completely, in 
the sense that they may even flow to two distinct non- 
compact phases (dense or dilute). In terms of the Li- 
ouville field theory of the FPL 2 model, one of the three 
bosons becomes massive and the other two decouple. 



2. Geometrical scaling exponents 

As discussed in Sec. || the complete spectrum of string 
dimensions may be computed within the interface rep- 
resentation from the two-point correlation functions be- 
tween height defects. Using the normalisation of Eq. ( fill ) 
the Burgers charge of the elementary vortices is now 
given by 

m 2 , = (2,0), m , 2 = (0,2), m x ,i = (1, 1) . (42) 

Furthermore, taking into account the compensating elec- 
tric charge needed to correct the spurious phase factors 
arising from the strings' winding around the defect cores, 
the electromagnetic charge of a general string defect gen- 
erating Sb black and s g grey strings becomes 

[e, m] Sb:Sg = [(7re b (l - <$ Sb ,o), 7re g (l - tf Sg , )) , (s b , %)]• 

(43) 

We note that there is now no distinction between an even 
and an odd number of strings. 

The string dimensions, obtained from Eq. (|3^), can 
finally be written as 



(e b , eg ) = < n < dl (eb)+< n ' dl (e g ) 



(44) 



where 



2x d - dl (e) = i(l T e) S 2 - T |^ 



(1 - 5 S , ) (45) 



is the formula for the string dimensions for one flavour 
of dense (dilute) O(rt)-loops. 



A. 0(n) model 



A particular variant of the O(n) model on the square 
lattice was defined by Blote and Nienhuis jlO) . Here the 
n-component spins live on the bonds of the square lat- 
tice whilst two and four-spin interactions are defined for 
spins that share a common vertex. This spin model has a 
graphical (loop) representation defined by the same ver- 
tices as the DPL 2 model. The only difference is in the 
assignment of vertex weights, and, more importantly, the 
grey bonds are simply treated as being empty. This 0(n) 
model is therefore a single-flavour loop model. We iden- 
tify the 0(n) loops with the black loops of the DPL 2 
model whose weight is = n. The weight of the DPL 2 
grey loops, which are formed by the empty bonds in the 
O(n) model, is n g — 1 (i.e., the weight of an 0(n) loop 
state does not depend on the number of loops formed by 
the empty bonds). 

If we focus for a moment on the central charge of the 
DPL 2 model along the line n g — 1, it is equal to the 
central charge of the black loops plus the central charge 
of the grey loops; see Eq. (41). The contribution of the 
grey loops to the central charge is or 1/2 depending on 
whether they are in the dense or the dilute phase; this, 
of course, is determined by the choice of vertex weights. 
As the black loops can also be either in the dense or di- 
lute phase this leads to four lines of critical points along 
which the central charge varies continuously. 
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Going back to the 0(n) model of Blote and Nien- 
huis, these authors identified, using numerical transfer 
matrix techniques, five lines of critical fixed points with 
exponents continuously varying with n, each specified 
by a particular choice of the vertex weights. One crit- 
ical line (branch 0) was directly mapped to the criti- 
cal (n + l) 2 -state Potts model, since along this line the 
flavour-crossing vertices 5 and 6 (see Fig. [l]) are excluded. 
Based on the numerical evidence (the measured central 
charge and a few of the scaling dimensions) two of the 
critical lines were identified with the usual dilute and 
dense phase of the O(n) loop model (branch 1 and 2), 
one as a superposition of the dense O(n) loop model and 
a critical Ising model (branch 4), whilst for the final fifth 
line (branch 3) the numerics were inconclusive. 

From the above analysis of the central charge it is clear 
that branches 1 and 2 can be identified with the grey 
loops being in the dense phase whilst the black loops are 
in the dilute and dense phase respectively. Along branch 
4 the black loops are dense whilst the gray loops are di- 
lute. This is also confirmed by comparing our predictions 
for the dimension xi.i (xh in Jl(J), Eqs. (|4|) and (45), 
with the numerically determined values in Ref. |hJ. We 
suspect that the "anomalous" branch 3 corresponds to 
both loop flavours being in the dilute phase. For n > 1 
this is born out by the numerics, but the correspondence 
seems to break down for smaller values of the black loop 
weight. Why this is so remains an interesting open ques- 
tion. 



B. Honeycomb lattice model 

We have mentioned earlier that the non-universality 
of the compact exponents forced us to base our theo- 
retical considerations on a particular lattice model. In 
order to examine the effects of having more than one 
flavour of loops we have opted for the square lattice. It 
is however reassuring to see that an altogether similar 
construction can be built starting from the honeycomb 
lattice, for which the fully-packed loop (FPL) model has 
been solved using Bethe ansatz methods S. 

The Coulomb gas construction for the honeycomb FPL 
model is based on the three-colouring of the lattice edges 
M . The colour vectors can be chosen as 




- ( ^A) , C- ( 



(46) 



and oriented loops are defined in the usual way as alter- 
nating sequences of colours B and C. In the correspond- 
ing densely-packed loop (DPL) model, where the vertex 
(A, A, A) is allowed and violates the three-colouring con- 
straint, the height field becomes one-dimensional: 



This is consistent with the magnetic plasma analogy, 
according to which the proliferation of defect charges 
m = 3A = (3, 0) would lead to the complete screening of 
the first height component in Eq. (46). 



ABC ACB 
BCA BAC 
CAB CBA 



ABC ACB 
BCA BAC 
CAB CBA 



ABC ACB 
BCA BAC 
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ABC 
BCA 
CAB 



0, 



B = 1, 



(47) 



1/3 1 4/3 2 7/3 3 

FIG. 7. Ideal state graph 1 of the DPL model on the hon- 
eycomb lattice. The repeat lattice is seen to be 1Z — Z. 

The FPL and the DPL models have the same six ideal 
states, and by projection onto the 2-direction in height 
space the DPL ideal state graph attains the appear- 
ance shown on Fig. |7]. Although each node now carries 
three labels the microscopic vertex weight A(x) remains a 
single- valued function onl. The repeat lattice is TZ = Z, 
and the electric charges live on the reciprocal lattice 
TZ* = 27rZ. For the background electric charge we find 
the unique solution eo = 7re, where n — 2cos(7re) is the 
loop fugacity. Since TZ* W — TZ* there are two candidates 
for the screening charges, effi — 2ir and effl = —2n. 
The former is associated with the most relevant vertex 
operator, and the loop ansatz fixes the (unique) elastic 
constant appearing in Se to be K = ^(1 — e). Finally, 
the central charge and the geometrical scaling dimensions 
are computed following the standard procedure, and as 
expected they are found to coincide with those of dense 
polymers. On the other hand, if a charge asymmetry 

(2) 

evolves e w is the correct screening charge, and the cor- 
responding solution for the elastic constant K = 5(l + e) 
is found to lead to the critical exponents of dilute poly- 
mers. 

It should now be clear that our ideas generalise to an 
arbitrary lattice of fixed coordination number z, provided 
that the usual Coulomb gas construction can be carried 
out. The description of the compact phase procedes via a 
z-colouring description, and the colour vectors C$ must 
satisfy Ylt=i Q = 0- Accordingly, the continuum field 
theory contains z — 1 massless bosons, which in the pre- 
sense of a background electric charge will couple in some 
non-trivial way. Now consider allowing violations of the 
fully-packing constraint. This will cause the [z/2] pairs 
of colour vectors, each defining one flavour of loops, to 
sum to zero separately, and, if z is odd, the remaining 
colour vector to vanish. The repeat lattice will then be a 
[z/2]-dimensional hypercube, and the loop flavours will 
decouple in the familiar way. Of the z—1 original bosons 
[(z — l)/2] become massive, whilst the remaining [z/2] 
describe each one independent flavour of dense or dilute 
polymers. 

A final remark pertains to our observation that in the 
DPL 2 model one of the two loop flavours may flow to the 
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dilute phase. Our Coulomb gas construction naturally 
suggests that it should also be possible for both flavours si- 
multaneously to be dilute. However, it is not obvious that 
this is compatible with the vertex configurations of Fig. [I], 
since the fractal dimension D[ = 2—X2 = (3+2e)/(2+2e) 
of dilute polymers is always less than two. 9 On the other 
hand, it is clear that the Coulomb gas construction of 



Sec. Ill would have remained exactly the same if we had 



allowed some additional DPL 2 vertices with two or four 
of the edges being unoccupied by any of the two loop 
flavours. Evidently, in this generalised model loops would 
be allowed to have any fractal dimension, 1 < Df < 2. 
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